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^ ■ Abstract 

Qh, As the first in a series of systematic work on dense hadronic matter, we study the properties of 

Q^l the pion in dense medium using Skyrme's effective Lagrangian as a unified theory of the hadronic 

interactions apphcable in the large limit. Dense baryonic matter is described as the ground 
state of a skyrmion matter which appears in two differentiated phases as a function of matter 
density: i) at high densities as a stable cubic-centered (CC) half-skyrmion crystal; ii) at low 
^ ■ densities as an unstable face-centered cubic (FCC) skyrmion crystal. We substitute the latter 

by a stable inhomogeneous phase of lumps of dense matter, which represents a naive Maxwell 
construction of the phase transition. This baryonic dense medium serves as a background for the 
pions whose effective in-medium Lagrangian we construct by allowing time-dependent quantum 
fluctuations on the classical dense matter field. We find that the same parameter which describes 
the phase transition for baryonic matter, the expectation value of the a field, also describes the 
phase transition for the dynamics of the in-medium pion. Thus, the structure of the baryonic 
ground state crucially determines the behavior of the pion in the medium. As matter density 
increases, {a) decreases, a phenomenon which we interpret to signal, in terms of the parameters 
of the effective pion Lagrangian /* and m* , the restoration of chiral symmetry at high density. 
Our calculation shows also the important role played by the higher powers in the density as 
it increases and chiral symmetry is being restored. This feature is likely to be generic at high 
density although our ground state may not be the true ground state. 



Pacs: 12.39-x, 13.60.Hb, 14.65-q, 14.70Dj 
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1 Introduction 



Understanding the properties of hadrons in a dense nuclear medium is currently an important 
issue in nuclear physics. The dynamically generated masses of the light hadrons reflect the 
symmetry breaking pattern of QCD. The quark condensate {qq) which represents the order 
parameter drops as the density of the nuclear medium increases. Consequently all light-quark 
hadron masses (except that of the pions) are expected to follow the behavior of the condensate, 
a basis of the scaling proposed in ^ and predicted in an effective field theory with the vector 
manifestation fixed point [2. The pion however is very special, because it is the Goldstone boson 
associated with the spontaneous symmetry breakdown of the chiral symmetry. The Hellmann- 
Feynman theorem 015) a niean field approximation in chiral perturbation theory [HI El 13 El j 
and the Nambu-Jona-Lasinio model [Sl llUlITT] lead to the same linear dependence in the baryon 
number density p namely, 

where S^^f is the pion-nucleon sigma term. If we assume that the Gell-Mann-Oakes-Renner 
(GMOR) relation holds for finite density, we have ^ 

f:^mf ^ -mg{qq)„ (2) 

where /* is the in-medium pion decay constant and m* the effective pion mass in the medium. 
There has been a lot of discussions on the physical consequence of eq. ^ with a decreasing 
quark condensate. The decrease in the pion decay constant is interpreted as a symptom of 
partial chiral restoration in the nuclear medium |10j. On the other hand, in refs. [HI El El) t^ie 
possibility of S-wave pion condensation resulting from the decrease of the pion effective mass as 
a consequence of the decrease of {qq) in eq. © has been discussed. However, if one lets the pion 
decay constant change, then the effective pion mass can even increase with increasing nuclear 
density as discussed in ref. [HI and a completely different scenario can then arise. 

The deeply bound pionic states observed in heavy nuclei (Pb, Sn etc.) at GSI ^1 El have 
renewed interest in these issues HH [inillZlllHl CHI EOlllIl ■ The well defined Is and 2p states 
with narrow width provide a new standard for the S-wave pion-nucleus optical potential Us{r), 
whose real part may be identified as the effective mass of the pions |14| II 5j 

mf{r) =ml + 2u;ReC/,(r). (3) 

This formula would seem to imply that the pion mass increases a 7% at the center of the Pb 
nucleus and a 3% in symmetric nuclear matter. But other possibilities exist. For instance, if 
one assumes the pion mass to be unmodified for low density which is not unplausible ^, then 
the measurements provide information on the scaling of the pion decay constant instead of that 
of the mass. See [T3l . 

Theoretical studies on the effective mass of the pion in a nuclear medium follow two different 
schemes. One traditional approach relies on optical potentials. However, it is difficult to relate 
the phenomenological vr-nucleus potential to the microscopic ir-N interactions in a systematic 



^For the moment, we are not distinguishing the time and space components of the pion decay constant. Which 
one we are dealing with wiU be specified later. 

^The pion mass of course remains strictly unmodified in the chiral limit but it is possible that even in the 
presence of small quark masses, the approximate chiral invariance protects the pion mass so that it is left more 
or less unaffected by density, at least at low density. An SU{2)c (2 color) QCD lattice calculation indicates that 
this is what happens |22|. 
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way. The other more modern approach proceeds via chiral perturbation theory [1^1 with exphcit 
baryon fields. By integrating out the baryon fields, using the mean field approximation, one 
obtains an effective theory whose parameters become dependent on the nuclear density p. These 
calculations consist of computing single-nucleon loops leading to a linear dependence on the 
density of the parameters. Higher orders in the nucleon loop expansion would provide us with 
higher powers of p, however, due to the strong couplings associated with the nuclear forces, it 
is not clear that the results will converge although chiral perturbation power counting is used. 
Perhaps more importantly an effective field theory Lagrangian properly matched to QCD at 
an appropriate scale as discussed by Harada and Yamawaki ^ must have a highly nontrivial 
"intrinsic" dependence on density (or temperature if in heat bath) in the parameters but this 
dependence is generally missing in the chiral Lagrangian so far used jSl I14j. 

We study the pion properties in dense baryonic matter in an approach based on Skyrme's 
Lagrangian, which we use to describe, in a unified manner, the dynamics of baryons, pions 
and the interactions between them. We proceed by calculating initially the ground state of 
infinite skyrmion matter. At high density it is described by a skyrmion crystal. At lower density 
the crystal suffers a phase transition and becomes unstable. The system chooses then as its 
ground state an inhomogeneous phase where skyrmions condense to form lumps of matter of 
higher density. We adopt a naive Maxwell construction to describe such state. We incorporate 
the pion in this medium by performing quantum fluctuations in the pion field of the Lagrangian 
describing this state. We obtain in this way the effective dynamics of the pion in dense skyrmion 
matter. The parameters of this Lagrangian are dependent on matter density and contain all 
higher orders in density as dictated by Skyrme's Lagrangian. 

Before proceeding we must clearly lay down the scope of this work here. We do not claim that 
the result we obtain - which is quite striking and intriguing - is in any way reflecting Nature. 
In fact it is perhaps likely that the "ground state" we obtain is not close to what it should be. 
More work is needed to obtain from what we have a true ground state and we have some ideas 
as to how to proceed. What we do here is that we will simply assume that we have certain states 
given by exact classical solutions of the field equations of a theory that is considered to be valid 
at large iVc, namely, the crystal structure of different symmetries, and study what they imply 
for the excitations thereon. This work should be taken as representing the first exploratory step 
towards a more realistic treatment - that is in accordance with QCD - of highly dense matter, 
the objective of our effort. 

This paper is organized as follows. In the next section, after a brief description of the model 
Lagrangian, we construct the classical skyrmion crystal solution using Fourier expansions. In 
sec. 3, we incorporate the pion fluctuations in the ground state and obtain the effective pion 
Lagrangian in the skyrmion medium. We discuss how the pion decay constant and the effective 
pion mass change in matter. Section 4 summarizes the most important results and contains our 
conclusions. 



2 The Model Lagrangian and the FCC Skyrmion Crystal 

The description of baryonic matter we are going to investigate is based on the Skyrme model 
Lagrangian [23j for the massive pion fields, which reads 

C = -i^Tt{U^df,UU^d''U) + ^Tip^df^U, U^d^Uf + :^^T¥(?7 + - 2), (4) 

where U G SU{2) is a nonlinear realization of the pion fields, /jr, e and m-,^ are the pion 
decay constant, Skyrme's parameter and the pion mass, respectively. Note that this Lagrangian 
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@ is described in terms of U , which contains only pionic degrees of freedom. We will not 
consider the isospin symmetry breaking due to the u- and d-quark mass difference, since it is 
of no relevance for the present study. Our aim is to study the properties of the pion in a dense 
baryonic medium. Since this is an effective Lagrangian, if matched to QCD j^j, its parameters 
should depend on the properties of the medium , i.e. matter density in our case as discussed in 
|24j . However we follow a rather different philosophy, namely we shall assume that the above 
Lagrangian embodies a complete description of all physical processes of interest, i.e., free pions, 
free baryons, many baryon states, dense baryonic matter and moreover, it also describes the 
pions interacting with any of these baryonic systems. We fix therefore the parameters of the 
Lagrangians in one of these processes, for example by describing the free pion and nucleon 
systems, and calculate how these values change as the density of the baryonic matter increases. 
From the recent understanding |21 123! effective theories that purport to match to QCD, we 
know that this is not totally consistent since EFTs can make sense only for given scales and the 
parameters of the theory necessarily flow. This issue which is quite intricate will be dealt with 
in a forthcoming publication. 

For completeness we should also mention that (jlj is the most economical effective Lagrangian 
that one can take. To be truly realistic, massive fields - such as the light-quark vector mesons, 
scalar mesons or "dilatons" ^ associated with trace anomaly etc - and/higher derivative terms 
need be incorporated but for our purpose it suffices to take the simplest yet sufficiently realistic 
form given by (0J). 

A similar procedure has been followed for example in chiral perturbation theory with explicit 
nucleon fields [71 |Hj . They start with a Lagrangian containing all the degrees of freedom and free 
parameters, they integrate out the nucleons in and out of an a priori assumed Fermi sea and in 
the process they get a Lagrangian density describing the pion in the medium which corresponds 
to the above Skyrme Lagrangian except that the quadratic (current algebra) and the mass terms 
pick up a density dependence of the form 

- T (^''^ + ^) '^(U^^^UU^d'^U) + (l - -p^p) Tr(t/ + U^- 2), (5) 

where p is the density of the nuclear matter and D^^^ and a are physical quantities obtained from 
the pion-nucleon interactions. Note that in this scheme, nuclear matter is assumed ab initio to 
be a Fermi sea devoid of the intrinsic dependence mentioned above. We shall discover in the 
next section that our free pion Lagrangian will be also modified in a similar manner but in the 
crudest approximation. At a more advanced level, there will be an intricate nontrivial density 
dependence missing in the chiral perturbative calculations. 

Having stated our scheme we proceed to describe dense baryonic matter in the Skyrme model. 
Consider a static field configuration Uq{x) = a+iT-TT with fj^+vf^ = 1. It is a nonlinear realization 
of the pion fields and a mapping, , from real space E? to the internal SU{2) space. If we take 
for Uq at infinity a constant SU(2) matrix, the configurations can be classified according to their 
winding number 

B = Jd^xpBix), (6) 

with 

PBix) = ^e^.kTriU^oViUoU^o'^jUoU^oVkUo). (7) 



^Implementing dilatons as was done in will play the role of incorporating the "intrinsic density dependence" 
inherent in Harada-Yamawaki theory of hidden local symmetry with the vector manifestation fixed point This 
point is being pursued and will be discussed in a forthcoming publication. 
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Skyrme |23j used this winding number to describe the baryon number. It is therefore through 
topology that this Lagrangian, written in terms of pion fields, describes baryonic systems. 

Various solutions of static field configurations which describe infinite matter can be obtained 
from skyrmions by assigning specific symmetries UHl IHl UHl H^l EOl • We proceed to describe 
in detail one of these solutions as a first step to study the behavior of the elementary pion 
in skyrmion matter. We adopt the method developed by Kugler and Shtrikman j^O]. They 
construct as an Ansatz for skyrmion matter a crystal whose configurations are such that each 
neighboring pair is rotated in isospin space, relative to the other with respect to the line joining 
them, by vr. Our starting crystal is FCC (face centered cubic) and the local baryon number 
takes its maximum at the FCC lattice cites. At high density, this configuration makes a phase 
transition to a CC (cubic centered) crystal, so called half- skyrmion" 27^ crystal, because only 
half of the baryon number carried by the single skyrmion is concentrated at the original FCC 
cites, while the other half is concentrated at the center of the links connecting those points, 
where the baryon number density was negligible initially. Thus the phase transition leads to a 
CC crystal made up of half-skyrmions. 

Let us construct the initial skyrmion matter solution. Consider a point in space x = (x, y, z) 
at which the field is given by (o", tti, 7r2, vra). Then, the FCC configuration is defined by the 
following symmetries: 

(51) under reflection in space x — > (— x, y, z), the field is also reflected in isospin space according 
to (fj, -vri,7r2,7r3), 

(52) under a rotation around the threefold axis in space x {y,z,x), the field is simultane- 
ously rotated in isospin space about the corresponding axis in isospin space according to 

{a, 7r2, TTs, TTi), 

(53) under a rotation around the fourfold axis in space x {x,z,—y), the field is rotated 
around the corresponding fourfold axis in isospin space according to (fi, 7ri,7r3, —7^2), 

(54) under a translation from a corner of a cube to the center of a face x — > {x + L,y + L, z), 
the field becomes rotated by vr about the axis perpendicular to the face according to 
[a, -TTi, -7r2,7r3). 

Here, 2L is the size of the single FCC unit cell which contains 4 skyrmions. Thus, the baryon 
number density is p = 1/2L^. The normal nuclear matter density po = 0.17/fm'^ corresponds to 
L ~ 1.43 fm. 

We initially introduce ^^unnormalized^ fields (a, vfi, 7f2, tts), which are then normalized^ to 
define properly C/q = a + ir • vf by 



with similar equations being satisfied by vrj {i = 1,2,3). The field configurations obeying the 
above symmetries can be easily found by expanding the unnormalized fields as Fourier series, 
i.e. 



a 



(8) 



a = 




a = 



Pabc cos(a7rx/L) cos{bTTy / L) cos{ctt z/L) 



(9) 



a,b,c 



and 



ahki sin(/i7rx/L) cos{kTTy / L) cos{Ittz/L) 

h,k,l 



(10) 
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(11) 



(12) 



The symmetry relations (^)-(^ll) restrict the modes appearing in eqs. (|9I12() as follows; 

(Ml) if h is even, then /c, I are restricted to odd numbers and a, 6, c are to even numbers, 

(M2) if h is odd, then fc, / are restricted to even numbers and a, 6, c are to odd numbers. 

Furthermore, the coefficients of the expansion satisfy othki = oihik and ^abc = l^bca = l^cab = 
Pacb = f^cba = Pbac- Note that the normalization process © does not spoil any symmetries that 
the unnormalized fields have, while the expansion coefficients Ohki and Pabc lose their meaning 
as Fourier coefficients in the normalized fields. 

Without loss of generality, we can locate the centers of the skyrmions at the corners of the 
cube and at the centers of the faces by letting a = —1 and 7rj(i = 1,2,3) =0 at those points. 
For the skyrmion field to have definite integer baryon number per cite, we should have a = +1 
and TTi{i = 1, 2, 3) = at the points such as (L, 0, 0). It gives the constraint 



If only the modes (M21) appear in the expansion, the configuration has an additional symme- 
try, namely that under the translation x (x + L,y, z), the field rotates under 0(4) by vr in the 
cr, vTi plane. This configuration corresponds to the half- skyrmion CC discussed above. Because 
of this additional symmetry, the physical quantities such as the local baryon number density 
and the local energy density become completely identical around the points with a = —1 and 
the points with a = +1. Thus, one half of the baryon number carried by a single skyrmion is 
concentrated at the cites where the centers of the single skyrmion is in the FCC crystal. The 
other half of the baryon number is now concentrated on the links connecting those points, where 
a takes the value +1, while in this location in the FCC configuration the baryon number density 
is low. As a consequence the average value {a) vanishes and we will see that this phenomenon 
signals that, in the dense medium, chiral symmetry restoration for the pion dynamics takes 
place. It is important to stress that it is the precise structure of the ground state which is 
responsible for the restoration. Both modes (Ml) and (M2) are included in the calculation. The 
half-skyrmion crystal arises as the stable ground state at high density, because the expansion 
coefficients associated with the modes (Ml) become suppressed. 

The Fourier series expansion method works well for the unnormalized fields and only a few 
modes are necessary.^ The expansion coefficients are used as variational parameters and de- 
termined by minimizing the energy of the configuration. The coefficients depend on the box 
size L. In Table 1, we list a few modes below E = 16(7r/L)^ and we use only the modes with 
E < 10(7r/L)^ in our calculation. 

In Fig. 1 we show the energy per baryon E/B as a function of the FCC box size parameter 
L. Each point denotes the minimum energy of the crystal configuration with a given value 
of L, which we find by using a proper minimization program such as "the down-hill simplex 
method" The solid circles correspond to the zero pion mass calculation and reproduce 

those of Kugler and Shtrikman [201 • In order for an easier comparison we present the results 

^Of course, through the normahzation process infinitely many modes appear in the normalized fields. 
However, the expansions converge rapidly. 




(13) 
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Table 1: A few modes used in the Fourier series expansion coefficients Ohki and Pabc in eqs. (PIT^ . 
E is the energy of the modes in units (vr/L)^ and d is the degeneracy of the mode. 
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with L given in units of (efT^)^^ (~ 0.45fm) and E/B in units of (67r^/^)/e (~ 1160 MeV), 
respectively. The latter enable us to compare the numerical results on E/B easily with its 
Bogolmoln'y bound for the skyrmion in the chiral limit, which can be expressed as E/B = 1 in 
this convention. 

By looking at the figure one sees that as we squeeze the system from L = 6 to around L = 3.8, 
the skyrmion system undergoes a phase transition from the FCC single skyrmion configuration 
to the CC half-skyrmion configuration. The system has a minimum energy at L ~ 2.4 with 
the energy per baryon E/B ~ 1.038, which is very close to the Bogolmoln'y bound. One can 
note in Fig. 1 that in the region L > 3.8, two different phases coexist. This coexistence leads to 
an interesting hysteresis phenomenon in the numerical process for finding the solutions, where 
one uses the found solution set for an L as the initial trial functions for the next L. If one 
starts from a small L, the numerical program might stay in a quasi-stable half-skyrmion CC 
configuration up to quite large values of L > 6, then suddenly it changes to the single skyrmion 
FCC configuration. This phase transition - and hysterisis - occurs in the unstable region where 
the system has negative pressure, so does not have any physical consequence. 

The phase to the left of the minimum is referred to in this work as ^^homogeneous" and there 
the background fields are described by a crystal configuration. The phase to the right of the 
minimum is ^Hnhomogeneous" since the pressure P = dE/dV is negative, the skyrmion matter 
is unstable, and to make it stable we will assume a realization in terms of an ^Hnhomogeneous 
phase" where the skyrmions condense to form dense lumps with large empty spaces. 

The open circles are the solutions found with a nonvanishing pion mass, m-rr = 140 MeV.^ 
Comparing to the skyrmion system for massless pions, the energy per baryon is somewhat higher. 
Furthermore, there is no first order phase transition at low densities. 

In Fig. 2, we present {a), i.e. the average value of a over space as a function of L. In the 
chiral limit, {a) rapidly drops as the system shrinks. It reaches zero at L ~ 3.8. where the 
system goes to the half-skyrmion phase. This phase transition can be interpreted, once the 
pion fluctuations are incorporated, as a signal for chiral symmetry restoration. In the case of 
massive pions, the transition in (a) is soft. Its value monotonically decreases and reaches zero 
asymptotically, as the density increases. Furthermore, as we can see in Fig. 3, where the local 

''Incorporating the pion mass into the problem introduces a new scale in the analysis and therefore we are 
forced to give specific values to the parameters of the chiral effective Lagrangian, the pion decay constant and the 
Skyrme parameter, a feature which we have avoided in the chiral limit. In order to proceed, we simply take their 
empirical values, that is, fw — 93 MeV and e — 4.75. Although the numerical results depend on these values, 
their qualitative behavior will not. 
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Figure 1: Energy per single skyrmion as a function of the size parameter L. The sohd circles 
show the results for massless pions and the open circles are those for massive pions. Note the 
rapid phase transition around L ~ 3.8 for massless pions. 



baryon number density is shown, for L = 2 (left) and L = 3.5 (right) in the z = plane, the 
system is (almost) a half-skyrmion crystal at high density. 



3 Pions in Skyrmion Matter 

Let Uq{x) = a + if ■ tt he the static skyrmion crystal configuration obtained so far. Then, we 
incorporate on it fluctuating time- dependent pion fields through the ansatz [SSj 

U{x,t) = VU^Uo{x)Vu^, (14) 

where 

= exp (^if ■ (j){x) / fj^ , (15) 

with (p describing the fluctuating pions. Substituting this ansatz into the Skyrme Lagrangian ^ 
and expanding the Lagrangian in terms of pions up to second order, we obtain the Lagrangian 
for the pion in the skyrmion matter which is given by 

C = CB{Uo) + CMiU^) + Ci, (16) 

where Cb{Uq) is the Lagrangian density for the skyrmion matter described by the static configu- 
ration Uq, Cm is for the Lagrangian density for the pions, and finally >C/ describes the interaction 
between the fluctuating pions and the skyrmion matter. Up to second order in the pion fields. 
Cm is simply the free pion Lagrangian to 

Cm = \d^.4>ad^<^a - \ml<g. (17) 
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£/ is quite complicated and is given by 



Hi = \4>aG''\x)4>b-\{d^(t>a)Hf^{x){dj4>b) 

+im2</,2Tr(l - Uo) + im2/,0„Tr(iT„C/o) (18) 



where 



C^^x) = ^Tr(^T^UoM-r'r'^+^^'Tr[[R^,r^][R^,A), (19) 
H^b^^S) = G^%, + ^^Tr(^[Ri,R,][f^,f']-[R„f'][R„r]^, (20) 

V^^ix) = ^f^Tr[iLi + Ri)Ta]+j^Tr(^[L,,Ta][U,Lj] + [Rj,Tam,R (21) 

A^{x) = ^/2Tr[(L,-i?i)r,] + ^Tr(^[L,-,rJ[L,,L,]-[i?,-,r„][i?i,i?,]) (22) 

with Li = {diU^)Uo and Ri = {diUo)Ul and f " = r'^ + [/or'^f/(|. The currents Vf" and 
produced by the static background field satisfy the following conservation equations: 

V-y" = 0, (23) 
= -\mlflTT{iTM. (24) 

The structure of our Lagrangian comes out to be very similar to that of chiral perturbation 
theory, eq. ©, and refs. [71 Let us use this fact to make some crude estimates. 

The interactions of pions and skyrmions as given by eq. (|18|) modifies the mass term in eq. (|17|) 

as 

\ml<j{x)<t^l. (25) 

This implies that the local effective pion mass can be even imaginary in the central part of the 
skyrmions where a is negative. By taking the average over space of o", we obtain an effective 
mass of the pion in the baryonic medium as 

ml ~ m^{a)^l\ (26) 

The kinetic term is also modified by the background, with respect to the free case, with an 
additional factor G°'^{x) given by 



G^'Hx) 



-{T^^Sab - VTaVTfe) + ^ { {{d^af + [S^^f^ {a^6ab + VTa^b) (27) 

diU-Ka - crdi-K^ {diCfT^b - (ydiTTf^ I . (28) 



The first term in G''*(f), -ir'^dab + vr aTTft, arises from the quadratic term of the Skyrme model 
Lagrangian, eq. Q, and corresponds to the term containing D^^ in eq. ©. Let us denote it by 
T"'^{x). It changes the kinetic term in eq. (|17|) to 

i/2(i + r"^(f))</.,<^fe. (29) 



9 



Ignoring that the inverse of this matrix T""^ becomes singular at the points where vr^ = 1 and 
taking the diagonal element, for example, F^^, eq. ()29() implies that the pion decay constant 
behaves locally as 

/:(x)~/.(i + rinf))^. (30) 



The average over the space leads to an effective pion decay constant in the medium /* given 
by 6 

/;~/.(l-|(vr2))i (31) 

If we absorb this change of 7,^ in the fields to keep the functional form of the theory as in the 
free case, eq. is modified to 

This equation shows the direct relation between the pion mass in the medium, m*, and the 
parameter (a) which signalled the phase transition in the case of skyrmion matter. We next 
show that chiral symmetry is restored in dense matter. 

In Fig. 4 we show the estimates of f^/ fn and m*/m7r as a function of the density. As the 
density increases, /* decreases to ~ 0.65/,r and then it remains constant at that value ^. Note 
that {a) 2 has the same slope at low densities, which leads to m^/m.,^ ~ 1 at low densities. Since 
at higher densities (vr^) becomes a constant, m*/m7r decreases like (o")^/^ with a factor which is 
greater than 1. As the density increases, higher order terms in p come to play important roles 
and m*/m7r decreases. 

The slope of {a) at low density is approximately 1/3. If we expand (a) about p = and 
compare it with eq. ©, we obtain 

(^)~l-^^ + ---l-|^P + ---, (33) 
3 Po f^rni 

which yields SttAt ~ m^/^/(3po) ~ 42 MeV, which is comparable with the experimental value 
45 MeV ^. This comparison is fully justified from the point of view of the jj expansion since 
both approaches should produce the same result to leading order in this expansion. The liner 
term is 0(1). 

The length scale is strongly dependent on our choice of the parameters f-^ and e. Thus one 
should be aware that the p scale in Fig. 4 could change quantitatively considerably if one chooses 
another parameter set, however the qualitative behavior will remain unchanged. 

We next proceed to a more rigorous derivation of these quantities using perturbation theory. 
The presence of the zero modes, associated with the collective motion of the skyrmion, and 
the complexity of G"'' make the quantization procedure non-trivial. Furthermore, G"^ becomes 
singular at the points where a{x) = 0. 

Let us obviate the complications associated with the zero modes and proceed. The momenta 
conjugate to the pion field cpa are given by 

U, = {dab + G'''m<Pb (34) 



^Comparing eqs. 1181 . 11911 and 1201 with eq. @ it is apparent that to second order Ci leads to diflerent time 
and space components for We keep the present discussion to leading order where they are equal. 

^These estimates are based on the lowest order approximation and should be valid only for low densities. 

*While this value is widely quoted, there is at present a considerable controversy on the precise value of this 
sigma term. In fact it can even be considerably higher than this. See |33| for a recent discussion. 
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0.00 



Figure 4: Naive estimates of ?7i*/m^ and f^/fn as functions of the baryon number density of 
skyrmion matter. 



which leads to the fohowing Hamiltonian 



c 



M 



(35) 



It is convenient, as in chiral perturbation theory j^, to decompose it into the free pion Hamil- 
tonian, TLq, plus the interacting Hamiltonian, Tii, defined by 



Ho = \Ma + \{^(t)a)-{^(pa) + \ml(t)a(t)a 

1 



1) 



(36) 



(37) 



Note that the momenta H^, defined by eq. (|34j) . and used in the derivation of the Hamiltonian 
do not appear in TCq. This simplifies matters considerably since only (j)a will define the canon- 
ical momenta of the quantization procedure. TCo defines the free pion propagator and TCi the 
interaction potentials as summarized in Fig. 5. 

We define the effective mass of the pion in skyrmion matter as the pole in the pion propagator 
in the limit of p — > 0, 

Gip) = 1 V (38) 

where S(|3) is the pion self-energy in matter. Thus, the interactions associated with the spatial 
derivatives acting on the pion fields will not play any role in the process up to second order in 
the interactions. Only the interactions characterized by G"'^ and a will be active to this order. 

Due to the crystal structure of the background fields in the homogeneous phase, G°'^{x) and 
a{x) can be expressed as 



<7 X 



E 



S{£) exp{i£-xTr/L), 



(39) 
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{Po,P) 
1 



(c) 



I? — 7)2 — 777,2 




Pq-P 



iPo,P) 



I (y: Sie)6'{p- q- 1) - s'ip- q) 

(d) 



{qo,q) 




ipo,p) 




ipo,p} 



-pl\EG^''i^)^'(p-^-^) 

(e) 




{po,p) 



-ieabcip + q)i\Yl VmS^ip- q- i) 



Figure 5: Free propagator and interactions between the pion fields and background skyrmion 
matter. 
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Table 2: Possible integer values for the expansion coefficients P^, 7^. and 7 i . 



/3f 


ix 


7II 


4 ^y 


4 


7!^ 


4 -^2/ -^2 


(i) 


even even even 


(i) 


odd odd 


odd 


(i) 


odd odd even 


(ii) 


odd odd odd 


(ii) 


even even 


even 


(ii) 


even even odd 



G''\x) = G"^(4exp(if-f7r/L), (40) 

with the help of a set of integers £ = (4i4'^^)- Thus, the corresponding a{p) and (^""^{p) in 
momentum space have the following structure, 

a{p) = J2 Si£) 5\p- i), r\p) = E G-\^¥\P- 4- (41) 

I i 

The FCC crystal yields non- vanishing expansion coefficients only for the integer set listed in 
Tab. 2. Furthermore, the coefficients have the following symmetries 

'S'(4; 4' 4) — ~l~'S'( 4; 4' 4) — ~l~'S'(4) 4' ^2;) — ~l~'S'(4) 4' 4) 
(4) 4' 4) — G ( 4) 4' 4) — (4; 4' 4) — ■I'G' (4) 4' 4) (42) 
(-^xj 4' 4) — G ( 4) 4' ^2) ~ (4) 4' 4) — ~I~G (4] 4' ^^) 

To the first order in the interactions, the self-energy gets contributions from cr(0) and G^^(O) 
given by 

S^(po,P - 0) = ml{S{Q) - 1) -p2Gii(0). (43) 
To this order, the effective pion mass m* becomes 



Since 5(0) is nothing but (fi), if we ignore the denominator, we can interpret our result as 
the naive estimate (j26p . Furthermore, keeping only the first order terms in G^^(x), the result 
coincides with our estimate ()32() . Thus, eq. 1)44(1 differs from eq. ()32() by the presence of the 
higher order term coming from the Skyrme term. Since the latter is positive definite, G^^(()) is 
always larger than — |(7r^) appearing in eq. ()32() . According to our numerical results, however, 
(_r^^(0) takes negative values for a wide range of p down to the density where E/B has the 
minimum. Thus, eq. (|44|) tells us that mJ^/m.,^ decreases slower than (c)^/^ in this range of 
but faster than the naive estimate H32|). 

If we include the second order diagrams as shown in Fig. 6, we obtain the self-energy as 

S^^(po,0) = Y.\po,Q) 

S{i)S{-i)mi 2S{i)G^\-i)mlpl 
^ \pI - ml- 47r2/L2 pl-ml- ^vr^/L^ 

(Gii(/)Gii(-/) + 2Gi^(4G^n-4 )vt \ 
^ pl-ml- 4^2/L2 J ■ ^ 



1/2 
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Figure 6: Diagrams used to evaluate the self-energy up to second order in the interactions. Here, 
b runs over 1,2,3 and ^ / 0. 

The resulting effective pion mass can be found numerically by solving the equation 

pl-ml-^"ipo,6) = 0, (46) 

after substituting pQ = m*^^ . 

In the chiral limit where = 0, eq. H46() always provides us with the trivial solution pQ = 0, 
which supports the consistency of our process. We show in Fig. 7 the numerical results on m'^/mn 
as a function of the density p. We have taken in eq. ()45() up to < 6 Fourier components. 
Although, there is no apparent reason for the interactions to be weak, m*^^ is very close to m*^. 
The Fourier expansion for S{p) and G"'^{p) can only be done for the homogeneous skyrmion 
crystal. 

Now that we have developed the perturbative expansion, we may check our estimate eq. (|31|) 
of the pion decay constant in the medium. We proceed by calculating it directly from the 
appropriate expectation value of the axial current. 

The total axial current of the system can be written as 

+l^eabcMx)y^^A^)+A^^Ax), (47) 

where A^^a and V^^a are the axial and vector currents from the background skyrmion matter, 
respectively. The total axial current satisfies the following conservation relation. 



d''A^,a = -Umla{x)<i)a{x) - ^€abcd^(l)b{x)V^A^) + ^-^mlTT{iTaUo). (48) 



The axial current of the fluctuating pion is obtained from eq. (|47j) by omitting the background 
(last) term of the equation. Let us use the symbol A'^^^ for the in medium pion axial current. 
Its conservation leads to the equation of motion, 

d^Af,^, = -Umla{x)Mx) - ^eabcd^4>b{x)V^A^). (49) 
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Figure 7: m*^ and m*J^ as a function of the density. We have used /7r=93MeV and e = 4.75 in 
the calculation. 



Applying the LSZ reduction formalism, the appropriate matrix element of the pion axial current 
in the ground state of skyrmion matter is, 

(0|<i(x)|vri(g) 0) = -i^ Uq,{l + G"(x))e-'?-, (50) 

where we have used vr^ to perform the calculation, exploiting isospin symmetry, and Z is the 
wave function renormalization constant obtained as = (1 + G^^(0)) in eq. (|4,3|) . The ground 
state of skyrmion matter is denoted by |0). In the spirit of our approximation we only consider 
the = component, to which does not contribute. If we regard the quantity appearing 
in the right hand side of the above equation as the in-medium pion decay constant and take the 
spatial average of G^^(x), we obtain, 

/; = /.(Vz)-i = ui + G^\^)Y'\ (51) 

which supports eq. (|31j) . Similarly, the matrix element for the conservation relation becomes 

(0|a^4i(:E)|7ri(g) 0) = -y/ZUmla{x). (52) 

Again, if we take the quantity in the right hand as —f*'^m^, this equation gives the following 
relation between free and in-medium quantities. 



<lJ1 = 

ml U 

Incorporating the relation eq. ()51|1. we obtain 

ml /2 



VZ(j(0). 



a(0). 



(53) 



(54) 
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Figure 8: (a) The homogeneous but unstable crystal phase and (b) the inhomogeneous but 
stable skyrmion matter. 



which is nothing but eq. (|32|) . 

Up to now, we have considered only the ^^homogeneous" phase where the background fields are 
in a crystal configuration. We next proceed to describe the ^Hnhomogeneous phase" . In Fig. 8, 
we present the situation schematically; in (a) the skyrmions are in the homogeneous crystal 
phase, while in (b) they condense to form rather dense lumps of matter and empty spaces. Our 
crude approximation consists of describing the dense lumps of matter by the crystal structure 
at the minimum energy per baryon ignoring surface effects. In this approximation E/B of the 
inhomogeneous phase (b) will be the same as {E / B)-cain- That is, the inhomogeneous matter 
has a lower energy than the homogeneous one and becomes stable. In Fig. 1, such energy is 
represented by a horizontal line starting at the minimum. Our approximation can be interpreted 
as a Maxwell construction between the two phases. 

In the inhomogeneous matter space is divided into two different phases. Let us distinguish 
these by using the subscripts "•yac" denoting the (empty) vacuum phase and ^^min" denoting 
the lumps of skyrmion matter at the minimum energy. The expectation values {a) can be 
approximated by 

, v(j) ^ (Q")miii^tnin + (q')vacKac /^^n 
''mm I 'vSLC 

where the superscript "h" and "i" denote that the quantities are for the homogeneous and 
inhomogeneous phase, respectively. 

In the inhomogeneous matter, the total volume of the dense droplets l^nin is nothing but 
L^jjj and the volume of the whole space is Vmm + K^ac- Thus, eq. can be simply written as 

(^)^^^ « (^)i^^ - - {atl){L,,,jLf. (56) 
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Figure 9: Maxwell construction for the homogeneous unstable crystal phase. We represent the 
change in (o") function of lattice size. 

The dash-dot lines in Fig. 9 represent the value of {a) for the inhomogeneous phase obtained 
by using the approximate Maxwell construction defined by eq. (|55|) . The change of {a) in the 
chiral limit is remarkable. Chiral restoration occurs at the point L = Lmin where E/B has 
its lowest value through a phase transition from the inhomogeneous phase to the homogeneous 
half-skyrmion CC phase at L = -Lmin- The sharp phase transition may be an artifact of our naive 
approximation (|55p . In the case of ^ also a rapid change in (a) occurs in the transition 
to the stable inhomogeneous phase, moreover there is almost no difference between the massive 
and massless cases. 

We show in Fig. 10 the value of m** obtained by eq. (jS2I) for the inhomogeneous phase. 

We learn from our calculation that the dynamics of the pion in the medium depends not 
only on the density of the baryonic matter but also from its detailed structure. It may be 
possible to interpret this as simulating the effect of the "intrinsic" density dependence required 
in Harada-Yamawaki theory |2l |2l] . 

4 Concluding Remarks 

We have studied the dynamics of the pion in a baryonic environment by incorporating quantum 
fluctuations to classical skyrmion matter. The procedure has led to an effective Lagrangian 
for the pion in the medium defined in terms of the parameters, /* and m*, which are density- 
dependent. The fact that we have a unified hadronic scheme to calculate both baryonic matter 
and the interactions of the pion in the background of such matter, with only a few parameters, 
allows us to draw important conclusions. This contrasts with the conventional approach where 
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the starting point is the matter-free vacuum, with the density effect taken into account pertur- 
batively. While our scheme is still far from being realistic, it contains the right symmetries and 
the support of physics, which renders it qualitatively predictive if not yet quantitatively. We 
are able to describe interesting phenomena though we are not able to predict yet the scales at 
which they will occur. 

The main result of our presentation is that, as the density of matter increases, chiral symmetry 
for the dynamics of the pion is restored, i.e. the effective parameters become those of a chiral 
symmetric phase. Two important dynamical mechanisms lead to this phenomenon. On the one 
hand the structure of the matter state is crucial in producing the phase transition. At high 
densities, the symmetries of the classical field describing the CC half skyrmion crystal, conspire 
to provide the a components with adequately distributed values on the lattice so that its average 
tends to vanish as the crystal becomes more and more precisely defined. This phenomenon 
which occurs at the level of the baryonic state is carried over by the quantum fluctuations to 
the dynamics of the pion in this baryonic medium. Thus the important physical quantities 
of the pion become density-dependent in terms of the expectation values of the classical fields 
and therefore are essential to this phase transition, which tends to restore chiral symmetry in 
the dynamics of the pion. On the other hand our procedure produces a density dependence 
which is non-trivial, i.e., the classical field leads to higher orders in the density dependence 
of the parameters. This intricate density dependence is analogous, if not entirely equivalent, 
to the intrinsic density dependence found in hidden local symmetry theory with the vector 
manifestation fixed point j24j . The results show that these additional terms, which do not appear 
in more conventional calculations, play an important role in the phase transition, suggesting that 
one should not trust results of this phenomenon obtained with the linear approximation or in 
few-order low-density expansions. 

It is important to stress that the results just exposed have required the development of a 
non-trivial perturbative procedure which has been solved beyond the leading order in the case 
of the effective pion mass. 

Ours is not yet a realistic calculation. The model we have used describes skyrmion matter, not 
nuclear matter. For instance we do not have the liquid structure of normal nuclear matter. In fact 
it appears nontrivial to arrive at a Fermi liquid from the crystal structure. Specifically while the 
in-medium pion decay constant going proportional to (o")^/^ up to near nuclear matter density is 
consistent with what one expects in the effective field theory with vector manifestation j24j . its 
saturation at higher density and non-vanishing as one approaches the chiral transition indicate 
that we may be describing a phase which is not the standard Wigner-Weyl symmetry. It may be 
in a pseudo-gap phase where the gap is non-zero though chiral symmetry is restored, resembling 
what might be happening in the "normal phase" of high Tc superconductivity ^- This may 
account for the fact that the in-medium pion mass decreases as one goes beyond normal matter 
density, which is at variance with the results of other approaches. 

It should be reiterated that the beauty of the whole procedure is that a unique Lagrangian 
is able to describe in a unified way all hadronic interactions, the ground state and fluctuations 
based on it. This hadronic theory provides us with results which are promising and exciting 
while exposing weaknesses of the more traditional approaches. Needless to say, there is still a 
long way to go before confronting Nature, such as, for instance, the properties of the pionic 
atoms that are offering a tantalizing view of what is going on in dense medium and of course 
the interior of compact stellar systems where chiral symmetry is presumably restored. Properly 
and realistically formulated, we might be able to "derive" from the flrst principles such novel 

'We are grateful to Kurt Langfeld for bringing our attention to this reference. 



19 



notions as BR scaling [Tl I24j . the hidden local symmetry/vector manifestation j2| etc. 
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